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Abstract 

Semidefinite  programs  are  a  class  of  optimization  problems  that  have  been  the  focus  of  in¬ 
tense  research  during  the  past  fifteen  years.  Semidefinite  programs  extend  linear  programs, 
and  both  are  defined  using  deterministic  data.  However,  uncertainty  is  naturally  present  in 
applications  leading  to  optimization  problems.  Stochastic  linear  programs  with  recourse  have 
been  studied  since  the  fifties  as  a  way  to  deal  with  uncertainty  in  data  defining  linear  pro¬ 
grams.  Recently,  the  authors  have  defined  an  analogous  extension  of  semidefinite  programs 
termed  stochastic  semidefinite  programs  with  recourse  to  deal  with  uncertainty  in  data  defining 
semidefinite  programs.  A  prominent  alternative  for  handling  uncertainty  in  data  defining  linear 
programs  is  chance- constrained  linear  programming.  In  this  paper  we  introduce  an  analogous 
extension  of  semidefinite  programs  termed  chance- constrained  semidefinite  programs  for  han¬ 
dling  uncertainty  in  data  defining  semidefinite  programs. 

AMS  Subject  Classifications:  90C15,  90C22,  90C05,  90C51 

Keywords:  Linear  programming,  stochastic  programming,  semidefinite  programming,  chance- 
constraints 


1  Introduction 

Semidefinite  programs  [1,  16,  18,  14]  are  a  class  of  optimization  problems  that  have  been  stud¬ 
ied  extensively  during  the  past  fifteen  years.  Semidefinite  programs  extend  linear  programs,  and 
since  both  are  defined  using  deterministic  data  we  shall  refer  to  them  as  deterministic  semidefinite 
programs  (DSDP’s)* 1  and  deterministic  linear  programs  (DLP’s)  respectively. 

Uncertainty  is  naturally  present  in  applications  leading  to  optimization  problems.  (Two-stage) 
stochastic  linear  programs  (with  recourse)  (SLP’s)  [6,  17,  7,  3,  8,  11]  have  been  studied  since  the 
fifties  as  a  way  to  deal  with  uncertainty  in  data  defining  DLP’s.  Indeed,  the  incorporation  of 
uncertainty  present  in  applications  into  models  is  so  important  that  stochastic  programming  is 
currently  one  of  the  most  active  research  subfields  of  optimization. 

‘Department  of  Mathematics,  Washington  State  University,  Pullman,  WA  99164-3113.  (ari@wsu.edu).  The  work 
of  this  author  was  supported  in  part  by  the  U.S.  Army  Research  Office  under  Grant  DAAD  19-00-1-0465  and  by 
Award  W911NF-08- 1-0530. 

^Division  of  Mathematical  and  Natural  Sciences,  Arizona  State  University,  Phoenix,  AZ  85069-7100. 
(yuntao.zhu@asu.edu).  The  work  of  this  author  was  supported  in  part  by  the  ASU  West  MGIA  Grant  2007. 

1  As  is  customary  in  mathematical  programming  and  optimization  literature,  we  use  the  term  deterministic  semidef¬ 
inite  program  to  mean  the  generic  form  of  a  problem,  and  the  term  deterministic  semidefinite  programming  to  mean 
the  field  of  activities  based  on  that  problem.  Both  will  be  denoted  by  the  acronym  DSDP,  while  the  acronym  DSDP’s 
will  denote  the  plural  of  the  first  usage.  Acronyms  DLP,  SLP,  SSDP,  CCLP  and  CCSDP  are  defined  and  used  in  the 
same  sense. 
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In  a  recent  paper  [2]  (see  also  [9]),  the  authors  have  defined  a  class  of  optimization  problems 
termed  (two-stage)  stochastic  semidefinite  programs  (with  recourse)  (SSDP’s)  to  deal  with  uncer¬ 
tainty  in  data  defining  DSDP’s.  SSDP’s  are  related  to  DSDP’s  in  the  same  way  that  SLP’s  are 
related  to  DLP’s. 

Chance- constrained  linear  programming  (CCLP)  [4,  5,  10,  11,  12]  is  a  prominent  alternative 
to  SLP  for  handling  uncertainty  in  data  defining  DLP’s.  In  this  paper,  we  introduce  a  class  of 
apparently  new  optimization  problems  that  we  refer  to  as  chance- constrained  semidefinite  programs 
(CCSDP’s).  CCSDP’s  are  related  to  DSDP’s  in  the  same  way  that  CCLP’s  are  related  to  DLP’s, 
and  CCSDP  is  an  alternative  to  SSDP  for  handling  uncertainty  in  data  defining  DSDP’s. 

The  rest  of  this  paper  is  structured  as  follows.  In  §2,  we  indicate  our  notation  and  some 
preliminaries  setting  the  stage  for  our  definition  of  a  CCSDP  given  in  §3.  The  main  reason  that 
semidefinite  programming  and  stochastic  programming  are  two  of  the  most  vibrant  research  fields 
of  optimization  is  their  applicability.  In  order  to  demonstrate  the  applicability  of  SSDP,  in  [2],  the 
authors  indicated  how  an  SSDP  may  be  formulated  for  dealing  with  uncertainty  in  data  defining 
the  minimum-volume  covering  ellipsoid  problem  [16,  13]  which  is  often  cited  as  an  example  DSDP. 
In  §4,  we  show  how  CCSDP  as  defined  in  §3  of  this  paper  provides  an  alternative  way  to  handle 
uncertainty  in  data  defining  the  minimum-volume  covering  ellipsoid  problem.  We  conclude  the 
paper  briefly  in  §5  commenting  on  how  the  definitions  in  this  paper  and  [2]  indicate  the  usefulness 
of  possible  joint  activities  by  researchers  in  semidefinite  programming  and  stochastic  programming. 


2  Preliminaries 


Our  notation  on  DSDP’s  follows  that  of  Todd  [14].  Let  Rmxn  and  RnVn  denote  the  vector  spaces 
of  real  m  x  n  matrices  and  real  symmetric  n  x  n  matrices  respectively.  For  U,  V  E  RnVn,  we  write 
U  0  (U  >-  0)  to  mean  that  U  is  positive  semidefinite  (positive  definite),  and  U  >z  V  or  V  A  U  to 
mean  that  U  —  V  >z  0.  For  U,  V  E  Kmxn,  we  write  U  •  V  :=  trace  ( UTV )  to  denote  the  Frobenius 
inner  product  between  U  and  V.  Given  C/j  E  RniXni  for  i  =  1, 2, . . . ,  n,  we  use  diag(C/i,  C/2, . . . ,  Un) 
to  denote  the  matrix  in  with  U\,  C/2, . . . ,  Un  on  the  diagonal  and  zeros  elsewhere. 

A  DSDP  in  primal  standard  form  is 

minimize  C  •  X 

subject  to  Ai  •  X  =  bi,  i  =  l,2,...,m  (1) 

x  y  0, 


where  A\  E  RnVn  for  *  =  1,2 ,m,  b  €  Rm  and  C  E  R”Vn  constitute  given  data,  and  X  E  R"Vn 
is  the  variable.  A  DSDP  in  dual  standard  form  is 

maximize  bT  y  ,  , 

subject  to  Ya= 1  ViAi  -  C 


where  A,  E  RnVn  for  i  =  1,  2, . . . ,  m,  b  E  Rm  and  C  E  RnVn  constitute  given  data,  and  y  E  Rm  is 
the  variable. 

It  is  possible  to  convert  a  problem  in  the  form  of  (2)  to  an  equivalent  problem  in  the  form  of 
(1)  and  vice  versa  (see  [18,  16]).  Also,  it  is  appropriate  to  refer  to  (1)  and  (2)  as  the  primal  form 
and  the  dual  form  respectively  due  to  the  duality  theory  that  exists  between  them  (see  [14,  §4]). 

A  DLP  in  primal  standard  form 


minimize  cTx 
subject  to  Ax 


> 


b 

0 


(3) 


2 


x 


and  its  dual 


maximize  b'y 
subject  to  AJy  <  c 

where  A  E  Mmxn,  b  E  Mm  and  c  E  constitute  given  data,  x  E  Mn  is  the  primal  variable, 
and  y  E  Mm  is  the  dual  variable  are  special  cases  of  (1)  and  (2)  respectively.  This  follows  from 
the  associations  C  :=  diag(ci,  C2, . . . ,  cn),  A{  :=  diag(Aji>  A$z,  ■  ■  ■ ,  Ajn)  for  *  =  1,  2, . . .  ,  m,  and 
X  :=  diag(xi,x2,  •  •  .,xn). 

Now  consider  the  DLP 

maximize  bT  y 

subject  to  AJy  <  c  (4) 

TJy  <  £ 

where  A  E  Mmxni,  T  E  M'mXTl2,  b  E  Mm,  c  E  Mni,  and  £  E  M™2  are  deterministic  data,  and  y  E  Mm 
is  the  decision  variable.  Clearly,  (4)  is  a  DLP  in  dual  standard  form. 

A  CCLP  based  on  (4)  is  defined  as  follows.  Suppose  that  the  vector  £  and  the  matrix  T  are 
both  random  depending  on  an  underlying  outcome  co  in  an  event  space  0  with  a  known  probability 
function  P.  Let  the  symbol  P  denote  probability.  Then  it  is  meaningful  to  require  that  the 
probability  of  the  constraints  TJ  {uo)y  <  £(cu)  being  satisfied  is  at  least  some  prescribed  value 
p  E  (0, 1),  rather  than  requiring  that  they  hold  for  all  outcomes  well,  This  leads  to  the  problem 

maximize  bTy 

subject  to  AJy  <  c  (5) 

^{TT(u;)y  <  £(w))  >  p 

which  is  termed  a  CCLP  [4,  5,  10,  11,  12]. 

Constraints  in  the  form  of  the  last  constraint  of  (5)  arise  naturally  in  various  applications  and 
are  called  chance-constraints  (or  probabilistic-constraints).  Such  constraints  can  be  viewed  as  a 
relaxation  of  the  requirement  that  constraints  are  enforced  for  all  possible  values  of  uncertain  data, 
which  could  be  prohibitive  or  even  impossible.  In  practice,  p  E  (0, 1)  may  be  chosen  close  to  1. 

3  Definition  of  a  CCSDP 

We  define  a  CCSDP  based  on  deterministic  data  A*  E  MniVni  for  i  =  1,2, ...  ,m,  b  E  Mm,  and 
C  E  MniVni;  and  random  data  Wt  E  Rn2Vn2  for  i  =  1, 2, . . . ,  m,  and  D  E  Mn2Vn2  whose  realizations 
depend  on  an  underlying  outcome  co  in  an  event  space  D  with  a  known  probability  function  P. 
Given  this  data,  we  define  a  CCSDP  as  the  problem 

maximize  bTy 

m 

subject  to  E  ViAi  A 

1=1 
m 

F{Y,yiWi(u)  p  d(u))  > 

i=  1 

where  y  E  Mm  is  the  variable,  and  p  E  (0, 1)  is  some  prescribed  value. 

4  An  application 

In  this  section  we  describe  an  application  and  its  formulation  as  a  CCSDP.  Our  application  is  a 
stochastic  version  of  the  minimum-volume  covering  ellipsoid  problem  (see  [16,  13])  often  cited  as 


C 

(6) 

p 
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an  example  DSDP.  In  [2]  we  formulated  an  alternative  stochastic  version  of  the  minimum-volume 
covering  ellipsoid  problem  as  an  SSDP. 


4.1  Preliminaries 

Suppose  that  we  are  given  ellipsoids  Si  C  Mn,  i  =  1,  2  defined  by 

Si  :=  {i  6  I"  :  xJ HiX  +  2 gj x  +  Vi  <  0}, 


where  Hi  £  MnVn,  Hi  y  0,  gi  £  Mn  and  Vi  £  M  for  i  =  1,  2.  Then  £i  contains  82  if  and  only  if  there 
is  a  r  >  0  such  that  matrix  inequality 


'  1 

T 

gi 

^  T 

'  h2 

T 

92 

L  9i 

n 

L  92 

V2  _ 

holds  [15].  Now  consider  the  ball  B  C  Mn  represented  by 

B  :=  {x  £  Mn  :  xTx  —  2xJ x  +  7  <  0}. 


The  center  of  B  is  x  and  its  radius  is  \J xTx  —  7.  The  distance  from  the  origin  to  the  center  of  the 
ball  is  V xTx. 

The  ball  B  contains  the  ellipsoids  S\  and  82  if  and  only  if  there  exist  t\  >  0  and  72  >  0  such 
that 

I  —x 
—xT  7 


<  Ti 

Hi 

T 

9i 

.  9i 

Vi 

i  =  1,2. 


4.2  The  Application 

We  describe  the  application  in  generic  terms  first.  Suppose  that  we  are  given  n/  fixed  ellipsoids 
Si  :=  {x  E  Mn  :  xJ HiX  +  2 gj x  +  Vi  <  0}  C  Mn,  i  =  1,2 ,...  ,nf.  Here  Hi  £  MnVn,  Hi  y  0,  gi  £  Mn 
and  Vi  £  M  for  i  =  1,2, ...  ,rif  are  deterministic  data.  We  are  also  given  nr  random  ellipsoids 
Si{u)  :=  {x  £  Mn  :  xJ Hi(u)x  +  2 gi(yo)J x  +  Vi(cu)  <  0},  i  =  1, 2, . . . ,  nr.  Here  for  i  =  1,  2, . . .  ,  nr, 
Hi(u>)  £  MnVri,  Hi((jj)  >~  0,  gi(io)  £  Vi{oj)  £  M  are  random  data  whose  realizations  depend  on  an 
underlying  outcome  to  in  an  event  space  Si  with  a  known  probability  function  P. 

Suppose  that  we  need  to  determine  a  ball  subject  to  two  types  of  constraints:  the  ball  must 
contain  all  n/  fixed  ellipsoids;  and  it  must  contain  the  nr  random  ellipsoids  with  probability  at 
least  p  £  (0, 1).  We  assume  that  the  cost  of  choosing  the  ball  has  two  components:  the  cost  of  the 
center  is  proportional  to  the  Euclidean  distance  to  the  center  from  the  origin,  and  the  cost  of  the 
radius  is  proportional  to  the  square  of  the  radius.  The  second  type  of  constraint  mentioned  above 
can  be  viewed  as  a  relaxation  of  the  requirement  that  the  ball  contains  all  the  realizations  of  the 
nr  random  ellipsoids,  which  could  be  prohibitive  or  even  impossible.  The  center  and  the  radius  are 
to  be  determined  so  that  the  total  cost  is  minimized. 

Before  proceeding  to  formulate  a  model  for  this  generic  application,  we  indicate  a  more  concrete 
version  of  it.  Let  n  :=  2.  The  fixed  ellipsoids  contain  targets  that  need  to  be  destroyed,  and  the 
random  ellipsoids  contain  targets  that  also  need  to  be  destroyed  but  are  moving.  Fighter  aircrafts 
take  off  from  the  origin  to  destroy  both  types  of  targets  with  a  planned  disk  of  coverage  having 
the  following  properties:  the  disk  contains  all  the  fixed  ellipsoids;  and  it  contains  the  realizations 
of  the  random  ellipsoids  with  probability  at  least  p  £  (0, 1).  Our  model  determines  the  center  and 
the  radius  of  the  disk  of  coverage  so  that  the  total  cost  is  minimized. 
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4.3  Formulation  of  the  Model 

Our  goal  is  to  determine  x  E  Mn  and  7  E  M  such  that  the  ball  B  defined  by 

B  :=  {x  E  :  xTx  —  2xJ  x  +  7  <  0} 


contains  the  fixed  ellipsoids  Si  for  i  =  1,2 and  the  realizations  of  the  random  ellipsoids 
Si(uS)  for  i  =  1,2 , . . .  ,nr  with  probability  at  least  p  E  (0, 1).  Note  that  the  center  of  B  is  x  and  that 
the  square  of  the  radius  of  B  is  xTx  —  7.  We  introduce  the  two  constraints 


d±I  x 
xT  d\ 


>-  0 


(7) 


and 


I  x 

xT  d2  +  7 


>-  0. 


(8) 


By  Schur  Complements  (7)  is  equivalent  to  d\  —  xT {d\I)~lx  >  0  which  in  turn  is  equivalent  to 
d\  >  V Constraint  (8)  is  similarly  equivalent  to  d2  +  7  —  xT/_1a;  >  0  and  to  the  constraint 
d2  >  xTx  —  7.  Thus  d]  is  an  upper  bound  on  the  distance  between  the  center  of  the  ball  B  and  the 
origin,  V xJx.  Similarly,  d2  is  an  upper  bound  on  square  of  the  radius  of  the  ball  B ,  xTx  —  7. 

Let  c  >  0  denote  the  cost  per  unit  of  the  Euclidean  distance  between  the  center  of  the  ball  B 
and  the  origin,  and  a  >  0  be  the  cost  per  unit  of  the  square  of  the  radius  of  B. 

We  define  the  decision  variable  x  E  ]gdTl/+Tl>’+» rH-3)  ag 


x  :  = 


d\i  d2,xJ , 


7,r 


where  except  for  the  auxiliary  variables  r  E  M nf  and  5  E  Mnr,  other  variables  are  as  specified  above. 
We  also  introduce  the  unit  cost  vector  c  E  R(n/+n’~+n+3)  as 


Then  we  get  the  model 

minimize 
subject  to 


c  :  = 


cTx 


c,  a,  0T,0,0T, 0T 


I 

-T 

—X 

^  Ti 

T 

9i 

—X 

7 

L  9i 

vi  _ 

,  i  1,2,...,  tij  , 


0  <  r, 

0  -< 


d\I  x 
xT  d\ 


(9) 


0  -c 


/  x 

xT  d2  +  7 


1  —x 
-xT  7 


^  di 
0  <5. 


Hi{io)  giiyj) 
,9i(v)  Vi(u)  _ 


,  i  =  l,2,...,nrSj  >p , 
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Problem  (9)  is  a  CCSDP  as  defined  in  §3.  We  now  indicate  the  assignments  that  need  to  be 
made  in  (6)  to  obtain  (9).  First  let  n\  :=  n/  +  nr  +  3 (n  +  1),  m  :=  n/  +  nr  +  re  +  3,  and  n2  :=  n  +  1. 
Then  let  y  :=  x  and  b  :=  —c.  We  use  0n  to  denote  the  zero  vector  in  MA,  and  0n  to  denote  the  zero 
matrix  in  MnVn.  Let  e*  be  the  z-th  elementary  vector  in  Mn,  and  In  be  the  identify  matrix  in  MnVn. 
Then  for  i  =  1, 2, . . . ,  (re/  +  nr  +  re  +  3),  we  define  matrix  Ai  G  ]R(n/+ri'-+3(n+i))v(n/+nr+3(n+i))  ag 
follows: 


-dl  diag(  0(n_)_i) ,  0 rifi  -^(n+l)>  ^(n+1))  0 nr  ); 


^2  •  diag^0(n+1),  Onjy, 

for  i  =  3, 4, . . . ,  (n  +  2) 

:=  diag 


On  On 

0 1  1 


j  0nr  )  J 


— e 


On 

T 


0-2) 


e(i—2) 

0 


)  Onj ) 


— e 


On 

T 


_e0— 2) 

0 


On  On 

0 1  1 


-4(n+3)  :=  diag 

for  i  =  (n  +  4),  (n  +  5), . . .  ,  (n/  +  n  +  3) 

4;  :=  diagf  -  [  9(i-n- 3) 


0-2) 
j  Onj)  0(n+l)) 


— e 


On 

T 


On  On 

o l  -i 


0-2) 

j  0nr  I  j 


—  e0“2) 

0 


,  diag(  e(j_n_3)),  0(n_|_3),  0(n_)_3),  0nr  J5 


9{i-n- 3)  v(i—n—3)  _ 

and  for  i  =  (re/  +  n  +  4),  (re/  +  n  +  5), . . . ,  (re/  +  nr  +  n  +  3) 

*d/  ■  diag(  0(n+i),  0(n_|_i),  0(n_|_3^,  diag(  e^_nj._n_ 3)))' 

Next  we  define  C  G  R(”-/+>v+3(n+i))v(n/+nr+3(n+i))  ag 


C  :=  diag 


In  On 

ol  0 


In  On 

ol  0 


,0, 


On^j  0(n+l)i 

Then  for  i  =  1, 2, . . . ,  (re/  +  nr  +  re  +  3),  we  define  matrix  Wi(u)  G  r0h-i)v0h-i)  as  follows: 

W\{lo)  :=  0(n+1); 

W2{u)  ■=  0(n+i); 

for  i  =  3, 4, . . . ,  (re  +  2) 


Wi(w)  :  = 


— e 


0-2) 


-e(Z-2)  0 


On  On 

ol  i 


^(n+3)M  :~ 

for  i  =  (n  +  4),  (re  +  5), . . .  ,  (rif  +  n  +  3) 

Wj(^)  0(n_|_3), 

and  for  z  =  (re/  +  re  +  4),  (re/  +  n  +  5), . . . ,  (re/  +  rer  +  re  +  3) 


Wi{u)  :=  - 


H(i— nf— n— 3)  0*0  9{i—nf—n— 3)(^0 

9(i-nf-n-3 )(W)  ^  (i-n  f -n-3)(^) 
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Finally,  we  define  D(oS)  G  M(n+1)v(n+1)  as 


D{u) 


In  0  n 

ol  o 


With  these  assignments  in  (6)  we  get  (9),  and  so  (9)  is  a  CCSDP  as  defined  in  §3. 


5  Concluding  Remarks 

In  this  paper,  we  have  defined  an  apparently  new  paradigm  for  stochastic  optimization  that  we  term 
chance-constrained  semidefinite  programming  (CCSDP).  CCSDP  is  an  alternative  to  (two-stage) 
stochastic  semidefinite  programming  (with  recourse)  (SSDP)  defined  in  [2]  for  handling  uncertainty 
in  data  defining  semidefinite  programs.  Our  developments  of  CCSDP  in  this  paper,  and  of  SSDP  in 

[2]  are  such  that  they  parallel  the  corresponding  developments  leading  to  chance-constrained  linear 
programs  (CCLP’s)  and  (two-stage)  stochastic  linear  programs  (with  recourse)  (SLP’s)  as  ways  of 
handling  uncertainty  in  data  defining  linear  programs.  As  a  consequence,  research  problems  on 
applications,  algorithms  and  theory  pertinent  to  CCSDP  and  SSDP,  and  possible  ways  for  their 
solution  essentially  suggest  themselves.  For  example,  it  readily  becomes  apparent  that  an  SSDP 
and  a  CCSDP  can  be  combined  into  a  single  model  similar  to  the  way  an  SLP  and  a  CCLP  is 
combined  in  [17]. 

Semidefinite  programming  and  stochastic  programming  are  two  of  the  most  vibrant  research 
subfields  of  optimization.  This  paper  and  [2]  suggest  exciting  new  possibilities  for  joint  work  by 
researchers  in  the  two  fields. 
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